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Abstract 

There has been mtense debate about the causes of the 1997-1998 
El Nino. One side sees the obvious intense westerly wind events as the 
main cause for the exceptional heating in summer 1997, the other em- 
phasizes slower oceanic processes. We present a quantitative analysis 
of all factors contributing to the onset of this El Nino. At six months' 
lead time the initial state contributes about 40% of the heating com- 
pared with an average year, and the wind about 50%. Compared with 
1996, these contributions are 30% and 90% respectively. As westerly 
wind events are difficult to predict, this limited the predictability of 
the onset of this El Nino. 



The Problem 

The 1997-1998 El Nino was one of the strongest on record. Unfortunately, 
its onset was not predicted as well as had been hoped (Pearce, 1997). In 
spite of claims that an El Nino could be predicted a year in advance, most 
predictions (Stockdale et al., 1998; Ji et al., 1996; Huang and Schneider, 1999; 
Kleeman et al., 1995) only started to indicate a weak event six months ahead 
of time. There have therefore been suggestions that El Nifio depends not 
only on internal factors, but also on external noise in the form of weather 
events in the western Pacific. 

The classical picture of El Nifio (Bjerknes, 1966; Philander, 1990) is that 
the usual temperature difference between the warm water near Indonesia and 
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the 'cold tongue' in the eastern equatorial Pacific causes an intensification of 
the trade winds. These keep the eastern region cool by drawing cold water 
to the surface. This positive feedback loop is kept in check by nonlinear 
effects. During an El Nino the loop is broken: a decreased temperature 
difference causes a slackening or reversal of the trade winds over large parts 
of the Pacific. This prevents cold water from reaching the surface, keeping 
the surface waters warm and sustaining the El Nino. 

This picture leaves open the question how an El Nino event is triggered 
and terminated. A variety of mechanisms has been proposed. On long time 
scales an unstable mode of the nonlinear coupled ocean-atmosphere system 
may be responsible (Neelin, 1991), either oscillatory or chaotic. Other au- 
thors stress the importance of a 'recharge' mechanism (Wyrtki, 1975; Jin, 
1997), with a built-up of warm water in the western Pacific preceding an El 
Nino. Another description on shorter time scales is in terms of refiections 
of equatorial Rossby and Kelvin waves in the thermocline (the interface be- 
tween warm surface water and the cold water below at about 100 m depth). 
These would provide the negative feedback that sustains oscillations (Suarez 
and Schopf, 1988; Battisti and Hirst, 1989; Kessler and McPhaden, 1995). 
However, short-scale atmospheric 'noise' in the form of westerly wind events 
in the western Pacific may also be essential in triggering an El Nino (Wyrtki, 
1985; Kessler et al, 1995). 

Here we trace the causes of the onset of last year's El Nino in May 1997 
over the six months from 1 December 1996. This is the time scale over 
which predictions are currently skillful. Although El Nino is an oscillation of 
the coupled ocean-atmosphere system, the analysis can be simplified by first 
studying the response of the ocean to forcing with observed wind stress and 
heat flux flelds. This response contains all time delays. The other part of the 
loop, the dependence of the wind stress and heat flux on the ocean surface 
temperature will be discussed separately. 

The ocean model used is the Hamburg Ocean Primitive Equation Model, 
HOPE (Prey et al., 1997; Wolff et al., 1997) version 2.3, which is very simflar 
to the ocean component of the European Centre for Medium-range Weather 
Forecasts (ecmwf) seasonal prediction system (Stockdale et al., 1998), but 
restricted to the Paciflc Ocean. It is a general circulation model with a hor- 
izontal resolution of 2.8°, increased to 0.5° along the equator, and a vertical 
resolution of 25 m in the upper ocean. It traces the evolution of temperature 
T, salinity S, horizontal velocities u, v and sea level 

This ocean model is forced with daily wind stress {Tx,Ty) and heat flux 
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Figure 1: The NlNo3 index observed (solid line) and simulated by the six- 
month forced model runs (dashed lines). 



Q from the ecmwf analysis, which in turn uses the excellent system of 
buoys (McPhaden et al., 1997) that observed this El Nino. Evaporation and 
precipitation are only implemented as a relaxation to climatological surface 
salinity. The initial state conditions are ecmwf analysed ocean states. To 
suppress systematic model errors we subtract a run starting from an average 
1 December ocean state forced with average wind and heat fluxes (both 
1979-1996 averages (Gibson et al, 1997)). 

The model simulates the onset of the 1997-1998 El Nifio quite well. We 
use the NINOS index A^3, which is a common measure of the strength of El 
Nifio (the anomalous sea surface temperature in the area 5°S-5°N, 90°W- 
150°W). In Fig. 1 the weekly observed NiNOS index (Reynolds and Smith, 
1994) is shown together with the index in the model run, compared to the 
same period one year earlier. The model overreacts somewhat to the forcing 
and simulates a NINOS index of 2.3 K at 1 June 1997, whereas in reality 
the index reached this value one month later. In 1995-1996 the simulation 
follows reality very well. 
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The Adjoint Model 

The value of the NlNo3 index at the end of a model run can be traced back 
to the model input (initial state, forcing) with an adjoint model. The normal 
ocean model is a (complicated) function M. that takes as input the state of 
the ocean at some time to (temperature Tq, salinity 5*0, etc.). Using the wind 
stress Ti and heat flux Qi for each day i for six months it then produces a final 
state temperature T„. The adjoint model (or backward derivative model) is 
the related function that takes as input derivatives to a scalar function of 
the final state, here the NIN03 index, dN^/dTn- It goes backward in time 
and uses the chain rule of differentiation (Gicring and Kaminski, 1998) to 
compute from these (and the forward trajectory) the derivatives ON-^/dTQ, 
dN^/dSo, dN^/dfi and ON^/dQi. These derivatives can be interpreted as 
sensitivity fields, giving the effect of a perturbation in the initial state or 
forcing fields. We can use them to make a Taylor expansion of the NIN03 
index to all the input variables: 




This means that the value of the index is explained as a sum of the influences 
of initial state temperature and salinity, and the wind and heat forcing dur- 
ing the six months of the run. These influences are each a dot product of the 
sensitivity to this variable (computed with the adjoint model) multiplied by 
its deviation from the normal state (extracted from the ecmwf analyses). 
To minimize higher order terms we take the average derivative from the sim- 
ulation and the climatology run. We have checked with actual perturbations 
that the accuracy of the linear approximation Eq. 1 is usually better than 
about 30% (within the model). Details can be found in van Oldenborgh et al. 
(1999). 

The 1997-1998 El Nino 

For the value of the NIN03 index on 1 June 1997 the linearization Eq. 1 gives a 
value of 1.8 K, compared with the 2.3 K simulated (and 1.3 K observed), this 
is within the expected error. The high value is mainly due to the influence 
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of the westerly wind anomalies (1.0 K) and the initial state temperature 
on 1 December 1996 (1.1 K). The salinity contributes —0.3 K, with a large 
uncertainty. 

The spatial structure of the influence of the initial state temperature is 
shown in Fig. 2. The top panel gives the temperature anomaly STq along the 
equator at the beginning of the run (Dec 1996), showing an unusually deep 
thermocline in the western Pacific and a shallower thermocline in the eastern 
Pacific. The second frame depicts the sensitivity of the June NIN03 index to 
temperature anomalies six months earlier, dN^/dTo. The third frame is just 
the product of the previous two; the integral of this over the whole ocean 
gives the 1.1 K contribution to the NINOS index mentioned before. The 
contribution is concentrated in the deeper layer of warm water along the 
equator in the western Pacific, in agreement with a 'recharge' mechanism. 

Fig. 3 shows the time structure of the influence of the zonal wind stress. 
The area under the solid graph gives the total influence, 1.0 K. The main 
causes of warming are the three peaks in zonal wind stress (dashed line) at 
the beginning of March, the end of March and the beginning of April, con- 
tributing about 0.6 K, 0.3 K and 0.5 K respectively. The peaks correspond 
with (very) strong westerly wind events in the western Paciflc. These gen- 
erated downwelling Kelvin waves in the thermocline that travelled east and 
deepened the layer of warm water in the eastern Paciflc 2-3 months later, 
increasing the surface temperature. There was also a strong wind event in 
December, contributing about 0.4 K over a negative baseline. From Fig. 3 it 
seems likely that it increased the strength of the later wind events by heating 
the eastern Paciflc in March. The heating effect of the March wind event also 
gave rise to an increase of the wind stress 5tx in May, but this reversal of the 
trade winds does not yet influence the NIN03 index dN^/dTx ■ St^, justifying 
the uncoupled analysis. 

The structure of the peaks in Fig. 3 can be seen more clearly in spatial 
views. In Fig. 4a the zonal wind stress anomaly 6tx is plotted for the second 
week of March. The westerly wind event corresponds to the large localized 
westerly anomaly around 150°E. Fig. 4b shows the sensitivity of the NiNOS 
index in June to the zonal wind stress during this week, dN^/dT^. This sen- 
sitivity consists of two main parts, both equatorally conflned. In the western 
and central Paciflc extra westerly wind stress would excite a downwelling 
Kelvin wave, raising the NiNOS index three months later. In the eastern Pa- 
ciflc the response would be in the form of a Rossby wave. The product of the 
anomaly and sensitivity flelds is shown in Fig. 4c. This gives the influence 



5 




dNs/dTo [K/K/sr/m] 




dNs/dTo ■ 5To [K/sr/m] 

Figure 2: Depth-longitude plots of the effect of the initial state temperature 
on the NlNo3 index in early June. At the top the analyzed temperature 
anomalies (averaged over 5°S-5°N) are shown at the beginning of December 
1996; the second frame depicts the sensitivity of the ocean to these tempera- 
ture anomalies and the third the product of these two, which gives the rise in 
the NlNo3 index on June 1 due to the thermal structure six months earlier. 
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Figure 3: The influence of the zonal wind stress on the NlNo3 index at 1 
June 1997 during the previous six months (sohd hne), the average anomalous 
wind stress over the area 130°E to 160°W, 5°S to 5°N. 

of zonal wind stress during this week on the NIN03 index, the integral of 
this field gives the corresponding value (0.22 K) in Fig. 3. The infiuence is 
contained in the intersection of the westerly wind event and the equatorial 
wave guide, and very localized in time and space. The effects of the other 
wind events are similar. 

The question remains whether the big influence of these wind events was 
due to their strength or to an increased sensitivity of the ocean dN^/dT,j.. 
We therefore repeated the analysis for the same months one year earlier, when 
the temperature in the eastern Pacific stayed below normal (Fig. 1). The 
adjoint model gives a NINOS index of —0.6 K, equal to the simulated index 
(the observed index was —0.7 K). This index is built up by a large negative 
infiuence of the wind stress, —1.5 K. and a positive influence of the heat flux, 
+0.9 K. The influence of the initial state temperature is also positive, but 
weaker than in the 1996-1997 +0.6 K, and the sahnity contributes —0.5 K. 

Although the built-up of warm water is also less pronounced, the largest 
difference is in the influence of the zonal wind stress. The sensitivity to 
zonal wind stress dN^/dr^ (over the area where its variability is largest) is 
compared for these two years in Fig. 5. During the time of the strong early 
March windevent the sensitivity was not very different bewteen the two years, 
but it was a factor two higher in April 1997 than in April 1996, and lower 
during the flrst two months. In all, these differences cannot explain more 
than a few tenths of a degree difference in the NIN03 index on 1 June. 

The difference between an El Nino in 1997 and no El Nino in 1996 can 
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Figure 4: The effect of tlie Marcli westerly windburst on tlie NlNo3 index 
in early June. At the top the averaged westerly wind stress anomaly for 
the week centered on 11 March 1997 is shown, the second frame depicts the 
sensitivity of the ocean to zonal wind stress and the third the product of 
these two which gives the rise in the NlNo3 index on June 1 due to this wind 
event. 
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Figure 5: The average sensitivity of the NlNo3 index on 1 June to westerly 
winds in the area defined in Fig. 3. 



be attributed for about 30% to an even stronger built-up of warm water in 
the western Pacific, and for about 90% to the the absence of strong westerly 
wind events in the western Pacific in the 1995-1996 rain season. A successful 
prediction scheme will have to predict the intensity of the westerly wind 
events correctly. However, the year-to-year variability of these wind events 
does not seem to depend on the state of the Pacific ocean (Shngo et al., 
1999), and at the moment is not predictable. 



Conclusions 

Using an adjoint ocean model we have shown that a successful prediction of 
the strong onset of the 1997-1998 El Nino, required a successful prediction of 
strong westerly wind events in March- April, which in our model contributed 
about 90% to the strength of the El Nino on 1 June 1997 compared to 
the situation one year earlier. The sensitivity to these wind events was not 
significantly different from the year before. The built-up of warm water 
contributed about 30% of the difference. The strong dependence on the 
westerly wind events would explain the relatively short lead time for correct 
predictions of the strong onset of this El Niiio. 
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